Fisher information manifestation of dynamical stability and transition to self-trapping 

for Bose-Einstein condensates 
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We investigate dynamical stability and self-trapping for Bose-Einstein condensates in a symmetric 
double well. The relation between the quantum Fisher information and the stability of the fixed 
point is studied. We find that the quantum Fisher information displays a sharp transition as the fixed 
point evolving from stable to unstable regime. Moreover, the transition from Josephson oscillation 
to self-trapping is accompanied by an abrupt change of the quantum Fisher information. 
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I. INTRODUCTION 

Quantum Fisher information (QFI) which character- 
izes the sensitivity of the state with respect to changes of 
the parameter, is a key concept in parameter estimation 
theory [H [2] . It has important apphcations in quantum 
technology such as quantum frequency standards [1, 0] , 
measurement of gravity accelerations [HI , and clock syn- 
chronization 0. Recently, the QFI was found to be able 
to detect entanglement and quantum phase transitions in 
many-body systems 0-[Il|- In the Lipkin-Meskhov-Glick 
model, the QFI can be used to characterize the ground 
state which displays a second-order quantum phase tran- 
sition [l^l. The critical point of a transverse Ising chain 
can also be estimated by QFI [ll[l3. 

We use QFI to study the dynamical stability and 
the transition from Josephson oscillation (JO) to self- 
trapping (ST) for Bose-Einstein condensates (BECs) in 
a double well potential. This transition is an interesting 
finding in BECs [l5l - l20j . By changing the atomic inter- 
action, the Josephson oscillation may be blocked, and 
the atoms of a BEC in a symmetric double-well poten- 
tial may show a highly asymmetric distribution between 
two wells. This phenomenon has been observed in the 
experiment pTj . 

Over the past few years, people found the transition 
is strongly related to the fixed points of the underly- 
ing classical dynamics. This is because quantum sys- 
tem can be mapped into a classical Hamiltonian which 
gives rise to fixed point solutions [l^]. Recently, some 
researchers have explored that quantum entanglement 
which is closely related to QFI can manifest the tran- 
sition [23r;25j . People also found quantum entanglement 
has relation to the classical fixed-point bifurcation— a loss 
of stability and the emergence of new fixed points (26l - 
[29| . However, effects of the stabilities of the fixed points 
on quantum information fiow is still lacking. It is thus 
certainly important to explore what features of the dy- 
namical stabilities are manifested by the QFI. 



We study the Hamiltonian of the system with the ap- 
proach of spin coherent state (SCS) associated to SU(2) 
group. We localize all the fixed points of the model in the 
phase space and discuss their stabilities. We also analyze 
the transition from JO to ST regime. In particular, we 
focus on the dynamics of the QFI for difi'erent interaction 
regions, such as ST and JO, stable and unstable regimes. 
We find that QFI can clearly demonstrate the dynamical 
stability and the transition from JO to ST. 

Our paper is organized as follows. In Sec. II, we discuss 
the parameter estimation, QFI, and introduce the maxi- 
mal mean QFI. In Sec. HI, we study the stabilities of the 
fixed points, and discuss the transition from JO to ST. 
Then in Sec. IV, we investigate the maximal mean QFI 
for two different initial states and reveal QFI manifesta- 
tion of the stability of the fixed point and the transition. 
Finally, a summary is provided in Sec. V. 



II. QUANTUM FISHER INFORMATION 

In this section, we discuss the QFI and the maximal 
mean QFI. Generally, for an input state pin under a lin- 
ear rotation by an angle Lp, the output state can be writ- 



ten as p^p — e 



Pine 



According to the Quantum 



Cramer-Rao theorem, the phase sensitivity ip has a lower 
bound limit [M [H 



1 



^JvF{pin,Jfi)'' 



(1) 



where ip is an unbiased estimator (i.e., (p — ^p), v \s the 
number of experiments, and F{pim Jn) denote the QFI 
which is defined as [30, 31] 



F{pin,Jn) = Tiip^L^j 



(2) 



Here, is the so-called symmetric logarithmic deriva- 
tive determined by the following equation 



(3) 
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In Eq. ([T]), one notices that, besides increasing the num- 
ber of experimental times v, we can improve the estima- 



2 



tion precision by choosing a proper state for a given 
Jfj, which maximize the value of the QFI. 

Now, we consider the maximal QFI for a given state. 
Based on Eq. ([3]), the expression of Eq. © is explicitly 
derived as 



(4) 

Pi + Pj 



where {\i)} are the eigenstates of with eigenvalues 
{pi} . Then the Eq. ^ can be compactly rewritten as 

F{pin,Jn)^riCn^, (5) 

where the matrix element for the symmetric matrix C is 

iP^-Pjf 



-'kl 



Pi + Pj 



mjk\j) {j\Ji\t) + mJi\j) (jVfcN)]. 

(6) 

It can be seen that the rotation along the ft direction 
affects the sensitivity of the state p. For a pure state, the 
QFI can be expressed as F{pin, Jn) = 4(A J^)^ [sl]- To 
obtain the maximal QFI, we rewritten the variance as 

{/\Jfif = nO{0^ CO)0'^ = n'Cdfi'^ , (7) 

where O is an orthogonal matrix, n' is a new direction 
defined as n' ~ nO, and is the diagonal form of C, 



Cd = O^CO==diag{Ai,A2,A3}, 



(8) 



where the A^'s are the eigenvalues of C. Now the maximal 
variance reads 

max(AJH)2 = max [\i{n[f + AzCr^)^ + hin'^f] . (9) 

In the above equation, the rotated direction is normalized 



and satisfies the condition 



= 1. If we set 



Amax — Ai as the maximal one of the eigenvalues, then 
il' = (1, 0, 0), and the original direction n = n'O^ . 
Now we get the maximal QFI 



For simplicity, we study the maximal mean QFI 



p — 

^ max — 



TV 



(10) 



(11) 



where N is the number of atoms. 



III. MODEL AND A CLASSICAL ANALOGUE 

System of BECs trapped in a symmetric double well 
have been well studied in theories and experiments (ssl - 
[36j . Due to the interaction between two wells, the system 
presents JO and nonlinear ST phenomena |18l - l20l |. For 
the two weakly coupled BECs system, the Hamiltonian 
can be described as i37h3al 



(12) 



where the angular momentum operators are defined in 
terms of the creation and annihilation boson operators 
a\ 2, ai,2 as 



(13) 
(14) 
(15) 



Jz 



a\a2 


+ 






2 




a{a2 




020.1 




2i 




a\ai 




02^12 




2 





which obey the SU(2) Lie algebra. The parameter il 
describes the coupling between two wells, and k denotes 
the effective interaction of atoms. In present work, we 
focus on the interaction strength k > and > 0. In 
the system, the total particle number N = a\ai + 0^20,2 
is a conserved quantity. 

To obtain the classical dynamics approach, we use a 
generalized SCS as an initial state, which is defined for- 
mally as [13411 



1^, 



-i9{Jx sin (jy — Jy cos <; 



3 , 



1/2 



m+j 



5 — I?, m) , (16) 



where j is the angular momentum quantum number, 

j = N/2, T = e""^tanf, and is the binomial 

coefficients. Under this SCS, the expectation values of 
the angular momenta are given by 

N 

{9 , (f)] J \6 , (j)) = — (sin 6 cos 4>, sin 9 sin 0, — cos 9) . (17) 

Meanwhile, we obtain the rescaled Hamiltonian (with 
constant terms dropped) 



H = {9,cf>\H\e,cb)/j 

= ft sin 6 cos 6 + Kr cos^ 



(18) 



where = {N — 1)k. With the help of path integral in 
the representation of SCS, we get the Lagrangian C for 
the system [isj 



C =N/2[h{l - cosO)^ - H], 



(19) 



which is associated with canonical coordinate and 
canonical momentum = H{1 — cos 6) For simplic- 
ity, we use p^ = —hcos9 as the canonical momentum. 
By setting h — I, the Hamiltonian can be rewritten 



I COS (f) + Kr 



(20) 



Then we obtain the following canonical Hamiltonian's 
equations of motion for p^ and in the phase space 



P4> = ~-P|)sin0, 
ilpcf, COS 4> 



— 2hirV(h 



(21) 
(22) 
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TABLE I: Fix points and stable regimes 



Parameters 


a 


b 


c 


d 




e = n/2, = 
stable 


9 = 7r/2, = TT 
stable 


N/A 


N/A 


Q < 2K,r 


e = Tv/2, = 
unstable 


e = 7r/2, = TT 
stable 


^ — arcsin[fi/(2/s:r)], = 
stable 


6 — n — arcsin[n/(2ft;r)], = 
stable 



It is noted that the motion governed by the above equa- 
tions are similar to that described by the mean field ap- 
proximation [isl [l9|. 

The stationary state solution of equations ((2T|) and 
can be obtained by assuming = (p — 0. For 
the case of = 0, we get = 1, (j) = 0, and = tt. For 

= 0, we obtain p^ = 0, p^ = y^l - cos 0)2, p^ = 

— y^l — ( 2^ cos 0)2". Under these conditions, we obtain 

several fixed points and list them in Table L It clearly 
shows that for stronger interaction, i.e. ft > 2Kr, in the 
phase space, there are two fixed points, corresponding to 
9 = tt/2, = and 6 — tt/2, cj) — tt, respectively. For 
weaker interaction < 2Kr, two more stable fixed points 

appear which correspond to = arcsin ^1 — ij^)^, 

= and 9 = TT - arcsin ^/l - (3777)^, = 0, respec- 
tively. 




1 2 0.5 1 1.5 2 

(|) /2ji <|) /27I 




0.5 1 1.5 2 0.5 1 1.5 



1^ /27I (|) /2jt 

FIG. 1: (Color online) Trajectories of the Hamiltonian system 
in the phase space for various parameters, corresponds to 
the double well population difference and represents the 
phase difference between the two wells. 



A. Dynamical stability 

In order to find out the fixed-point bifurcations, we 
need to analyze the stabilities of the fixed points. We 
first discuss the two fixed points {9 = 7r/2, = 0,7r). 
These points are interesting because they depend on 
neither tunneling strength nor self-collision interaction 
strength. We shall adopt the linear stability analysis that 
has wide applications in variable nonlinear systems. To 
begin with, we assume P4, — + 5p^ and = 0" + J^, 
where (p^ , 0" ) denote one of the fixed points in the phase 
space, 5p^ and represent the deviations in the popula- 
tion difference and relative phase from the fixed points, 
respectively. With Eqs. (PT|) and (1^^ . we obtain the lin- 
earized equation 



dt\d. 



M 



where M is the Jacobian matrix 



M 



dpsdct> 




(23) 



(24) 



The eigenvalues of the linearized equation may be real, 
pure or complex imaginary. As is well known, the eigen- 
values of the Jacobian matrix depict the types and sta- 
bilities of the fixed points. By substituting 9 = 7r/2, 



= 0, vr into Ec^. ([M)) and calculating the eigenvalues, we 
can find that the fixed point 9 = 7r/2, = is stable for 
Q, > 2Kr- When ft < 2Kr, it becomes unstable. While for 
the fixed point 9 — 7r/2, = tt, it is always stable. Com- 
pany to the condition that two new fixed points emerge 
for ri < 2k J., the classical bifurcation condition then can 
be obtained as 



(25) 



The stabilities of the fixed points can also be seen from 
the trajectories of the Hamiltonian. In Fig.(IT]), we plot 
the evolution trajectories in the plane of 9 and for var- 
ious tunneling strengths. It shows that a fixed point is 
stable if the evolution trajectories are loops around a 
fixed point; otherwise it is unstable. Note that all the 
stabilities of the fixed points which obtained by the nu- 
merical simulation are consistent with the above theoret- 
ical analysis. 



B. Dynamical transition between JO and ST 
regimes 

Now we consider the transition from JO to ST regime. 
In the JO regime, the population difference oscillates 
symmetrically between two wells, and its average is zero. 
In the ST regime, the average of the population differ- 
ence is nonzero, which can be obtained by the following 
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FIG. 2: (Color online) Fidelity of the overlap for different 
interaction of the two-mode BEC. The blue curve represents 
a weaker interaction with Q = Kr, whereas the red curve with 
dots represents a stronger interaction with Q = 4k,. . For both 
cases the initial state of the system is the SCS \d ~ = O) 
and the number of atom is A'' = 100. 



condition [l^ 



il sin ^0 cos (po + Kr cos^ 9o > 



IV. EFFECT OF STABILITY ON QUANTUM 
DYNAMIC 

First we investigate the effect of stability on quantum 
dynamic of the system. We start with the initial state 
g,s 9 — 7r/2, (f) — 0, (i.e. which corresponds to 

a fixed point. Such a state can be realized by applying 
a two photon tt/2 pulse to the state \j,j) with all the 
atoms in the internal state |F = 1, nip — 1) [46l - [50| . For 
this initial state, in classical analogue, it will change its 
stability at the classical bifurcation. 



A. Effect of stability on Fidelity 

In treating the quantum dynamical problem, it is help- 
ful to bear in mind some results from quantum infor- 
mation theory concerning fidelity. Fidelity has been 
widely studied for the problems of transition probabil- 
ity in quantum mechanics. Mathematically, for two pure 
states, the fidelity is defined as 



F = |(*| $)| . 

For the system with the initial wave function 
state at arbitrary time t can be expanded as [5! 




(26) 



=I]c™(t) |j, m), 



(29) 



with 9o and (po being the initial condition. From Eq. 
([26l) . we get the critical point of the transition 



1 — sin 6*0 cos (f>o 



(27) 



For the initial value — 0, (po — 0, the transition 
parameter is flc — Kr which is consistent with the re- 
sult of the model obtained by the mean field approxi- 



mation [H, [Hi 

4>o = TT and 6*0 



For 9o = 7r/6, </>o = 0, fie = ^. If 



7r/6, the corresponding critical value 
becomes flc = It shows that the critical point can be 
adjusted by the relative phase between two wells, and the 
relative phase can be experimentally adjusted by using a 
"phase-imprinting" method [45| . 

We emphasize that the particle number is large in the 
above classical analysis, however, in the practical exper- 
iment, the number of particles is finite. In the following 
discussions, we consider the quantum dynamic of the QFI 
to investigate the dynamical stability of the fixed point 
and the transition from JO to ST. We will consider two 
cases that the initial SCSs are chosen as 16* = 7r/2, = 0) 
and |6l = 0, (/) = 0). The SCS \9^0,<j) = 0) is simply a 
Dicke state \j, —j) in which all atoms lie in one of the 
wells. \9 — tt/2,(P ~ 0) is a state with a well-defined 
phase difference </> = and is a phase state of a two- 
mode boson system. 



and the amplitudes Cm(i) obeys 



icmit) = 2Km^Cm(t) -|- —Cm~i{t)y/{j + m){j - m + 1) 

n 



+ -Cm+i{t)^{j -m){j+m + l). (30) 



With the help of Eq. ([16)) , the amplitude for the initial 
SCS can be calculated as 

c™(o)=^(„?+,y^'. (31) 

According to the Heisenberg function, we get 



Jz 



(32) 
(33) 
(34) 



Unfortunately, for this model, it cannot be solved ex- 
actly for a many-particle case. Numerical results of the 
wave function overlap between the initial and the evolved 
state F = K7r/2,0| ^'(<))| are plotted in Fig. dSJ. The red 
curve with dots correspond to £7 = and blue one cor- 
respond to ~ Kr- For n — AKr, it can be seen that the 
fidelity oscillates around the value 1 . It indicates that the 
dynamic does not take the state far away from the initial 
state. According to the stability analysis in Sec. Ill, the 
fixed point which corresponds to the initial state is stable 
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FIG. 3: (Color online) Top panel shows interaction-time plot 
of maximal mean QFI Fmax- The bottom shows the -Fmax 
as a function of Q/itr with rescaled time (a) nt — 6 and (b) 
Kt — 24. Initial state of the system is SCS \6 — ^,<j> = O'^ and 
the number of atom is A'^ = 500. 



in the regime > 2kj.. When il = k^, we can find that 
the fidelity oscillates between and 1, and the amplitude 
of the oscillations is inhomogeneous. In the classical ana- 
logue, the fixed point is unstable in the regime fl < 2^^. 
From the above analysis, we find that there is a perfect 
classical-quantum correspondence. 



B. Effect of stability on QFI 

Now, we consider the effect of stability on the QFI. 
Numerical results of the maximal mean QFI i^max as a 
function of and Kt are plotted in Fig. (jSj. Obviously, 
the dynamic of the quantum system can be well illus- 
trated by the maximal mean QFI. As is shown in the top 
panel of Fig. ^ , the behaviors of the maximal mean QFI 
is quite different on the two sides of the classical bifurca- 
tion point n — 2k,-. To clearly depict the phenomenon, 
at the bottom of Fig. we plot Fmax as a function of 
fl/Kr for two arbitrary rescaled time nt = 6 and Kt — 24, 
respectively. These two figures show that Fmax behaves 
irregular oscillations in the regime < 2Kr, in which 
the fixed point is unstable. When > 2Kr, the fixed 
point is stable and we can find Fmax oscillates with time 
around the initial value. Moreover, we can find that the 
value of the maximal mean QFI in the stable regime is 
much smaller than that in the unstable regime. The phe- 
nomenon can be understand by the fidelity which have 
been discussed in the above. In the stable regime, the 
evolved state is close to the initial SCS, and the F^ax is 
small. While in the unstable regime, the evolved state far 
away from the initial SCS, and the -Fmax become bigger. 
These results suggest that the QFI can be well used to 
characterize the stability of the fixed point in this model. 




0.2 0.4 0.6 0.2 0.4 0.6 



Kt Kt 

FIG. 4: (Color online) Quantum mechanical time evolution 
of (Jz) in the ST and JO regimes for the initial condition 
\j, —j), corresponding to all atoms initially in one of the wells. 
The plots in terms of a dimensional parameter Kt, with (a) 
Q = 3Kr; (b) n = 2Kr/3. The number of atom is = 100. 



V. DYNAMICAL TRANSITION BETWEEN JO 
AND ST REGIMES 

A. population difference in different regimes 

Below, we consider another dynamical problem, the 
transition from JO to ST regime. Firstly, we investigate 
the population difference between two wells. We choose 
\9 — 0,(j) — 0) as the initial state, which corresponds to all 
atoms localized in one well. In this case, it is clearly that 
for f2 = ^Kr the initial state is related to the JO regime, 
whereas for = 2Kr/3 it is related to the ST regime. In 
Fig. (|4^) and Fig. dHa), we plot the quantum evolution of 
(Jz) as a function of Kt for the JO and the ST regimes, 
respectively. In the JO regime, {Jz) oscillates around 
zero during the evolution. It indicates that there is no 
preferential tunneling to any of the wells. While in the ST 
regime, {Jz) oscillates around a non-zero value and only 
in the half plane. It shows that part of the condensate 
is trapped in one of the wells. From the above results, 
we find that the dynamical properties of such a quantum 
system are quite different for the JO and the ST regimes. 



B. QFI in different regimes 

To well understand the quantum dynamical transition 
from JO to ST phenomenon, we calculate the time evo- 
lution of the maximal mean QFI. For a pure Dicke state 
I J, m) , the maximal QFI is obtained as 

F„,ax = 2(j2+j-m2). (35) 

Then for the initial state \j, — j), the maximal mean QFI 
can be got as flnax = 1- Numerical results of the maximal 
mean QFI -Fmax with different interactions are plotted 
in Fig. dS]). From the top of Fig. dS]), it is clearly seen 
that the behavior of the maximal mean QFI are quite 
different in the JO and the ST regimes. To well describe 
the behaviours of the maximal mean QFI for different 
regimes, we plot the maximal mean QFI as a function of 
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FIG. 5: (Color online) Top panel shows interaction-time plot 
of maximal mean QFI Fmax- The bottom shows the -Pmax 
as a function of Q./hir with rescaled time (a) nt = & and (b) 
K,t = 25. Initial state of the system is the SCS |6 = 0, </> = 0) 
and the number of atom is iV = 500. 

fl/Kr with two rescaled time Kt — 6 and Kt = 25 in the 
bottom of Fig. ([S]). From these two figures, we can see 
that, for < Kr, the maximal mean QFI increases slowly. 
As the interaction strength approach the transition point 
ilc = Hr, the -Fmax incrcascs quickly. When > k^, it 
oscillates quickly and its value is much larger than that 
in JO regime. It is shown that the QFI can be served 
as an indicator of the transition from JO to ST regime. 
Such abrupt change of the QFI also can be found in the 
phase transition of spin systems . 

In experiments, so far, many observers have found 
number squeezing in this model (53| . Moreover, some 
researchers have shown that the number fluctuation has 



a nontrivial relation with the phase fluctuation, even for 
a single- mode light held |5^]. We hope that the transi- 
tion of the maximal mean QFI for this system can be 
observed in experiment in the future. 



VI. CONCLUSION 

In conclusion, we have studied the dynamical evolu- 
tion of a two mode BEG in a symmetric double well. 
All the relevant flxed points and their stabilities were 
analyzed. We investigated the dynamical behavior of 
maximal mean QFI in this system which well described 
the stability of the flxed point. For the initial state 
\6 = n/2,(f) ~ 0), the numerical result showed that -Fmax 
has a small oscillation around the initial value in the 
51 > 2Kr regime, while in the regime fl < 2Kr, we found 
that -Fmax displays a strong irregular oscillation. 

We also investigated the maximal mean QFI in the 
JO and the ST regimes which corresponds to the differ- 
ent distribution of particles in two wells. For the initial 
state 1^ = 0, = 0), we found the mean maximal QFI in- 
creases quickly at the critical point which corresponds to 
the boundary between the JO and the ST regimes. From 
these results, it showed that the QFI not only character- 
izes the stability of the fixed point, but can also signal 
the presence the transition from JO to ST regime. 
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